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Abstract 

The modular vector field plays an important role in the theory of Poisson manifolds and is intimately 
connected with the Poisson cohomology of the space. In this paper we investigate its significance in the 
theory of integrable systems. We illustrate in detail the case of the Toda lattice both in Flaschka and 
natural coordinates. 
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1 Introduction 

The infinitesimal generator of the modular automorphism group is important in the cohomology 
theory of Poisson manifolds. On the other hand, one of the main characteristics of integrable 
Hamiltonian systems is the formation of hierarchies consisting of Poisson tensors. Each one of 
these tensors carries naturally a modular vector field. One may, therefore, speak of a hierarchy 
of modular vector fields. The purpose of this paper is to begin an investigation of the role of 
this sequence of vector fields in the theory of integrable systems. As an illustrative example, 
we choose the famous Toda lattice, a system that is one of the most basic in the theory of finite 
dimensional integrable models and closely related to the theory of simple Lie groups. There is 
no doubt, however, that the results of this paper hold for other similar systems as well. 

The modular vector field appears first in a paper of Koszul [17], as a special case of an 
operator on contravariant differential forms of degree —1. The same vector field was later used 
by Dufur and Haraki in [8] to classify quadratic Poisson brackets in IR 3 . It was called "curl", 
which in K 3 is dual to "divergence". In fact, the operator of Koszul applied to a vector field, 
gives the usual divergence of differential calculus with respect to the standard volume form. 
One can therefore speak of the divergence of a Poisson bracket; this is what the modular vector 
field is. The divergence of a Poisson bracket is useful in the classification of Poisson structures 
in low dimensions, e.g., [19| [T3j. 

One class of operators considered by Koszul is the following. Let A p be the space of covariant 
antisymmetric tensors on a smooth orientable manifold M of dimension n and let A p denote 
the space of covariant differential forms. Choose a volume form uj and define an operator 

D : Ap > Ap-i 

which satisfies 

D 2 = 0. 
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The volume form uj induces an isomorphism $ from A p to A n p . Define the modular operator 
as: 

D = r 1 odo$, 

where d is the exterior derivative. 

Weinstein gives a different interpretation. Suppose one has a Poisson tensor tt and a 
smooth positive volume form uj. Consider the operator X u , which sends a smooth function / 
to div^Xj, where Xf is the Hamiltonian vector field generated by / with respect to tt. It turns 
out that this operator is a derivation and hence a vector field. It coincides with the Koszul 
operator, D, at the level of contravariant 2-tensors. It satisfies: 

L Xbj u = 0, L Xuj tt = 0. (1) 

Even though we deal with spaces that are orientable, we remark that in the non-orientable 
case one can either replace volumes by densities or use the recent approach of [14). 

In local coordinates ( ), the modular vector field of the Poisson tensor tt is given 

by the formula: 

w = >:o:^i^-- (2) 

Lichnerowicz [18] considers the following cohomology defined on the space of contravariant 
tensors of a Poisson manifold. Let (M, tt) be a Poisson manifold. Define a coboundary operator 
d w which assigns to each p-tensor A, a (p + l)-tensor d n A given by 

d w A = -[tt,A], 

where [•, •] denotes the Schouten bracket. We have that d^A = [tt, [tt, A}} = and consequently 
defines a cohomology. More details on the Schouten bracket can be found in [TH [221 127] . 
An element A is a p-cocycle if [tt, A] = 0. An element B is a p-coboundary if B = [tt, C], for 
some (p — l)-tensor C. Let 

Z n (M, tt) = {A : [tt,A] = 0} 




and 

The quotient 



B n (M, tt) = {B : B = [tt,C]} 

Z n (M,TT) 



H n (M, 



TT 



B n (M, tt) 

is the nth cohomology group. The elements of the first cohomology group are the infinitesimal 
automorphisms of the Poisson structure modulo the Hamiltonian vector fields. It follows from 
(JT|) that the modular vector field is an element of the first cohomology group. If one replaces uj 
by au), where a is a positive smooth function, then 



Xau) — X w + X\ na ■ (3) 

Therefore, the modular vector field is a well-defined element of the first cohomology group in 
the Lichnerowicz cohomology. Since it forms an element of the cohomology group, Weinstein in 
[29] uses the term "modular class". It makes sense to consider two such vector fields as equal 
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if in fact they differ by a Hamiltonian vector field. The reader can also refer to [TS1 [2E1 EZJ for 
more details on Poisson manifolds and cohomology. 

The modular operator D is a graded derivation. The following relations hold for a general 
2-tensor n and vector fields X, Y: 

D[n,X] = [D(n),X]-[7v,D(X)] (4) 
D[X,Y]=XD(Y)-YD(X). (5) 

We note the minus sign in (j3J) that appears because we use the convention XJ(g) = {/, g}m in 
the definition of the Hamiltonian vector field XJ . 

The purpose of this paper is to present a complete study of the modular hierarchy in the 
case of the Toda lattice. In Section [2] we review the finite, classical Toda lattice and its multi- 
Hamiltonian nature. In Section [3] we examine the modular hierarchy in Flaschka coordinates 
(a, b) G M 2Ar_1 . In particular, we establish the Hamiltonian character of the modular vector fields 
associated to the well-known hierarchy of Poisson tensors for the Toda lattice, see Theorem [3J 
and present some of their basic properties in Corollary [TJ A new bi-Hamiltonian formulation of 
the Toda lattice is provided in Theorem [2j In Section H] we study the infinite modular sequence 
in natural coordinates (q,p) G M. 2N , and present a formula that iteratively produces all members 
of this sequence. In Section [5] we comment on the results of this paper and present them in 
compact form. 

2 The Toda lattice 

The Hamiltonian of the Toda lattice is given by 

N N-l 

H( qi ,...,q N , Pl ,...,p N ) = J2 ^Pl + Y, e " +1 ■ ( 6 ) 

i=i i=i 

This type of Hamiltonian was first considered by Morikazu Toda [26]. Equation ([6]) is known 
as the classical, finite, non-periodic Toda lattice to distinguish the system from the many and 
various other versions, e.g., the relativistic, quantum, infinite, periodic etc. The integrability of 
the system was established in 1974 independently by Flaschka [IT], Henon [TJ] and Manakov 
|21j. The original Toda lattice can be viewed as a discrete version of the Korteweg-de Vries 
equation. It is called a lattice, as in atomic lattice, since interatomic interaction was studied. 
This system also appears in cosmology and the work of Seiberg and Witten on supersymmetric 
Yang-Mills theories. It has applications in analog computing and numerical computation of 
eigenvalues. However, the Toda lattice is mainly a theoretical mathematical model which is 
important due to the rich mathematical structure encoded in it. 
Hamilton's equations take the form 

Qj = Pi 

p j = fflj-i-u - e ii-<ij+i ; j = i ) . . . j n. 

The system is integrable. One can find a set of independent functions {Hi, . . . , H^} which are 
constants of motion for Hamilton's equations. To determine the constants of motion, one uses 



3 



Flaschka's transformation: 



i e §(9i-9i+l) 



h = —zPi, 



1 = 1, 



N-l 
N. 



The equations of motion become 



a-i (k+i - h) 

2(«?-ii)- 



(7) 



(8) 



These equations can be written as a Lax pair L = [B, L], where L is the Jacobi matrix 



and 
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This is an example of an isospectral deformation; the entries of L vary over time but the 
eigenvalues remain constant, i.e. Aj = 0. It follows that the functions Hj = jtrL^ are constants 
of motion. We note that 

N 

H 1 = X 1 + X 2 + --- + X N = J2 b i 

i=l 

corresponds to the total momentum and 

N N-l 



i=l 



1=1 



is the Hamiltonian. 

Consider R 2N with coordinates (91, . . . , qN,Pi, ■ ■ ■ ,Pn), the standard symplectic bracket 



N 



If I =Vf— — -— — 
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and the mapping F : M. 2N — > M. 2N 1 denned by 

F : (q u . . .,q N ,px, . . . ,p N ) -> (a 1: . . . ,a/v-iA, • • • , M- 

The Flaschka transformation F is a symplectic realization of a degenerate Lie Poisson 
bracket on IR 2 ^" 1 , i.e. there exists a Poisson bracket on M 27V_1 which satisfies 

{f,g}oF = {foF,goF} s . 

This bracket (up to a constant multiple) is given by 

{cii,bi} = -at 
\Q>i, bi+i} = di ; 

all other brackets are zero. We denote this bracket by 7Ti. Its Lie algebraic interpretation can 
be found in [16]. In iri, the only casimir is Hi = b\ + b 2 + ■ ■ ■ + b^, and the Hamiltonian is 
H2 = \ tr L 2 . The invariants Hi are in involution with respect to 7Ti. For a proof of these facts 
see [I]. 

The quadratic Toda bracket appears in conjunction with isospectral deformations of Ja- 
cobi matrices. Let A be an eigenvalue of L with normalized eigenvector v. Using standard 
perturbation theory one obtains, 

VA = (2v l v 2 , 2v N ^v N , v 2 ,..., v 2 N Y := U x , 
where VA denotes (J^-, • • • , d ^ 1 , ■ ■ ■ , Some manipulations show that U x satisfies 

n 2 U x = \tv 1 U x , 

where n\ and 712 are skew-symmetric matrices. The defining relations for the Poisson tensor 112 
are the following quadratic functions of the a« and bi : 

{ a ii a i+l} = 2~ a i a i+l 

{di,bi} = -dibi 
{bi,b i+1 } = 2d 2 ; 

all other brackets are zero. The quadratic Toda bracket appeared in a paper of Adler [TJ in 
1979. It is a Poisson bracket in which the Hamiltonian vector field generated by Hi is the same 
as the Hamiltonian vector field generated by H2 with respect to the iri bracket. The tensor 112 
has det L as casimir and Hi = tr L as the Hamiltonian. The eigenvalues of L (and therefore 
the Hi as well) are still in involution. Furthermore, 1^2 is compatible with tti. We have 

n 2 VH i = >K X VH j + x , J = 1,2,... . (10) 

These relations are similar to the Lenard relations for the KdV equation; they are generally 
called the Lenard relations. Taking j — 1 in ffTU]) . we conclude that the Toda lattice is bi- 
Hamiltonian. Bi-Hamiltonian structures were introduced by Magri in [20] • Using results from 
[3], ones proves that the Toda lattice is multi-Hamiltonian: 

n 2 VHi = 7Ti V# 2 = TTO V# 3 = 7T_! V H 4 = ... . (11) 
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The Hamiltonian hierarchies of the Toda lattice are well-known. The results are usually 
presented either in the natural (q,p) coordinates or in the more convenient Flaschka coordinates 
(a,b). In the former case the hierarchy of higher invariants are generated by the use of a 
recursion operator [7J [10] . We remark that recursion operators were first introduced by Olver 
[25] . The system is bi-Hamiltonian and one of the brackets is symplectic. Thus, one can find a 
recursion operator by inverting the symplectic tensor. The recursion operator is then applied 
to the initial symplectic bracket to produce an infinite sequence of Poisson tensors. However, in 
the case of the Toda lattice in Flaschka variables (a, b), the first two Poisson brackets Hi and 7r2 
are non-invert ible and therefore this method fails. The absence of a recursion operator for the 
finite Toda lattice is also mentioned in Morosi and Tondo [23], where a Ninjenhuis tensor for 
the infinite Toda lattice is calculated. The family of Poisson tensors in this case is constructed 
using master symmetries. Invariant functions and Hamiltonian vector fields are preserved by 
master symmetries. New Poisson brackets are generated using Lie derivatives in the direction 
of these vector fields, and they satisfy interesting deformation relations. We quote the results 
from references [3J, H] . 

Theorem 1 There exists a sequence of vector fields X i; for i > —1, and a sequence of con- 
travariant 2-tensors ttj, j > 1, satisfying : 

i) 7Tj are all Poisson. 

ii) The functions Hi, % > 1, are in involution with respect to all of the nj. 
hi) XiiHj) = (z + j)H i+j , i > -1, j > 1. 

iv) Lx.iTj = (j - i - 2)Tv i+j , i > -I, j > 1. 

v) [X t , X j ] = (j-i)X i+j ,i>0, j>0. 

vi) 7Tj V Hi = Kj-i V H i+ i, where ttj denotes the Poisson matrix of the tensor Ttj. 
Remark 1: Theorem 3 was extended for all integer values of the index in [5]. 

3 The modular class in Flaschka coordinates 

We consider the modular vector field of the Poisson tensor nj, denoted by Yj = D(jTj), where 
D is the Koszul operator. We begin with some preliminary results needed to prove the main 
theorem that establishes the Hamiltonian character of the modular class. The vector fields of 
Theorem [U denoted by Xj, are master symmetries for the Toda lattice system. For example, 
X\ is given as follows: 




N-l 



N 



where 



A = -ia-ih + {i + 2)aib i+1 

B % = (2i + 3)a 4 2 + (1 - 2i)a?_! + b}. 



If we define the function / = ln(ai • • • ajy_i), then we have the following proposition: 
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Proposition 1 

D(X 1 )=X 1 (f) + 2H 1 . 

Proof. 

By definition 

- E§f + Ef 

i=l i=l 

n— 1 n 

= X)H&i + (* + 2)6^1]+ 2 

i=l i=l 

= X 1 (/) + 2Tr(L) = X 1 (/) + 2/f 1 . 
A similar relation holds for the second master symmetry X 2 that is defined as 

N- 1 a N a 

*-E<+E<, 

1=1 1=1 



where 



d = (2 - i)a?_!ai + (1 - z)ai&- + aibib i+1 



i-l 



+ (i+ l)aia? +1 + (i + l)aj6- +1 + a? + (J^ 6 J -)a i (6 i+ i - &;) 

i-l i-2 



3=1 3=1 
+ (3 - 2i)at 1 b l -i + (4 - 2z)a 2 _ 1 + 6?. 

Proposition 2 

D(X 2 ) =X 2 (f) + 6H 2 . 

Proof. 

Similar to the proof of Proposition [TJ 

We can generalize the results of Propositions [T] and [2] as follows: 
Proposition 3 

D(X n ) =X n (f)+n(n+l)H n . 

Proof. 

The proof is inductive. The result holds for n = 1,2. We assume that the result holds for n > 2 
and prove it for n + 1 . 
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D(X n+1 ) 

n — i 

(X 1 D(X n ) - X n D(X 1 )) 

(X x (X n (f) + n(n + l)H n ) - X n (X x (f) + 2H X )) 

([X x ,X n ] + n(n + l)X x (H n ) - 2X x {H n )) 

((n - l)X n+x (f) + (n(n + 1) - 2)X x {H n )) 

= X n+X (f) + (n + 2)X x {H n ) = X n+X (f) + (n + \){n + 2)H n+x . 

Below we prove a statement that we will use in Theorem [31 but it is also important in its 
own right. It is a new bi-Hamiltonian formulation of the Toda lattice. 

Theorem 2 

ATjg = Xp+\ where g = ln(det(L)), and j > 1. 

Proof. 

We use the Lenard relations for the eigenvalues [HE], ttj V A, = A, ; 7Tj_iV Aj. 

^Hx - ^-(Ai+.-.+Ajv) 



ttj-VCAi H h A 



A' J 



On the other hand 



- ^lnAi...Ajv 



7Tj +1 V In Ai ... A 



N 



N 1 



i=l 

The following theorem investigates the Hamiltonian character of the divergence of the infi- 
nite sequence of Poisson tensors 7tj. In particular, it states that the divergence of the Poisson 
tensor ttj is the Hamiltonian vector field given by the function h := ln(ai • ■ ■ a/v-i) + (j — 
l)ln(det(L)). 

Theorem 3 

Yj = Xj+y_ x \ g , where f = ln(ai • • -ajv-i), and g = ln(det(L)). 



8 



Proof. 

We will prove the proposition inductively in two steps. First we will show that it holds for 
j = 1 and then for 2 < j < 3. Consequently, we will show that it holds for j : = 4 and then for 

i>5. 

We have that Y\ = D(tti), where n\ is given by 

{ai,bi} = -a { 
{di,b i+1 } = -di. 

Thus, using (T5]) we obtain 

Y ^w 1 -wr {0 '-- ' 1 ' -' -- 1) '- 

On the other hand, it is not hard to show that the Hamiltonian vector field induced by / = 
ln(ai • • -a^v-i) with respect to the bracket tti has the form: 

a, = 0, i = 1, . . . , N — 1 

b\ = 1, b N = -1 

bi = 0, i = 2, ...,N- 1. 

Hence Y x = Xf . 

Now, let us assume that the theorem holds for j > 4. We will then show that it holds for 
j + 1. By definition L Xl ^j = [Xl,7Tj] = [tt^Xi] and using {iv) of Theorem [Q we have that 
^i+i = (Fsyh'^i]- Thus, 

Y j+1 = D(7r j+1 ) = -^D([n j ,X 1 }) 

([D(7T j ),X 1 ]-fa,D(X 1 )]) 

([Y j ,X 1 ]-fa,D(X 1 )}) 

{[xj\x 1 ] + (j - i)[x;\x 1 ] - fa, d(x 1 )]) 

{[[Ltt^X,] + (j - l)^ 1 ,^] - fa,D{X 1 )]) 

(fa, [X u /]] + [/, fa, X,}] + (j - l)(j - 2)*g - [vr,, D(X x )]) 

(fa,X 1 (f)} + (j - 3)[/, 7T i+1 ] + (j - l)(j - 2)*g - [TTy, D(Xi)]) 

([vr,, D(X X )] - 2[7r i; Fx] + (j - 3)^ +1 + (j - l)(j - 2)Ag - [vr,, D(X 1 ))) 
(j-l)(j-2)-2 
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In an identical manner one proves that the theorem holds for 2 < j < 3. To show that the 
proposition holds for j — 4 we use a similar argument as the one used above, however, we 
employ the relation 7r 4 = — \L X ^2 = — |[^2)7r2], to 

Corollary 1 

a) Y i (H j ) = Y j (H i ) 

b) LyJXj = —Ly^i- 

Proof. 

For part a) we have: 

Y t {H 3 ) = {X? + (i- 

= 

= {f, Hj} Vi 

= Vfir.VH, 

= VfcVHi 

= {f,H i } 7Tj 

= Yj(Hi). 

For part b), 

[[h, Tlj], 7Tj] + [[TTj, -Ki\,h] + [[iTi, h],TTj] = 0. 

Since [vrj,^] = the result follows. 

4 The modular hierarchy in (q,p)— coordinates 

The bi-Hamiltonian structure of the Toda lattice appears in a paper of Das and Okubo in 1989 
[7]. The generation of master symmetries and Poisson tensors using Oevel's theorem and the 
connection with the results in Flaschka coordinates is due to Fernandes |1Q| . In principle, the 
method is general and works for other finite dimensional systems as well. For example, this 
approach was used by Nunes da Costa and Marie [2] in the case of the relativistic Toda lattice. 
The procedure goes as follows. One defines a second Poisson bracket in the space of canonical 
variables (qi, . . . , qN,Pi, ■ ■ ■ ,Pn)- This gives rise to a recursion operator, call it 71. The presence 
of a conformal symmetry, as defined by Oevel, allows one to use the recursion operator and 
generate an infinite sequence of master symmetries. These, in turn, project to the space of 
the new variables (a, b) to produce a sequence of master symmetries in the reduced space. We 
quote the relevant theorem of Oevel that appears in [24]. 

Theorem 4 Suppose that Zq is a conformal symmetry for both Poisson tensors J\, J2 and 
function hi, i.e. for some scalars X, fi, and v we have 

£-z J\ = AJi, £z J2 = A^2, £>z hi = vh\. 
Then the vector fields Zi = 7Z l Z are master symmetries and we have, 

(a) C Zl h k = (u + (k - 1 + - X))h i+k 
(6) C Zi J k = + (k - % - 2)(/i - X))J i+k 
{c)[Zi, Z k ] =(/jt — X)(k- i)Z i+k . 



avoid division by zero. 



W 1 ) W 
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We note that master symmetries were introduced in [12] . 

We proceed with the computation of the modular vector fields that are associated to the 
family of Poisson tensors given in (q,p)— coordinates. The first Poisson tensor in the hierarchy 
is the standard canonical symplectic tensor: 




where ijv denotes the N x N identity matrix. The second Poisson tensor has the form, 

J = ( ^ N Bn\ 
2 \— B N Cn J 

where An is the N x N skew-symmetric matrix defined by = 1 = — a^, for % < j, = 
diag(— pi, —p2, ■ ■ ■ , ~Pn), and Cat is the NxN skew-symmetric matrix whose non-zero terms are 
given by c^+i = — = e qi ~ qi+1 , for % — 1, 2, . . . , N — 1. If we let Ji = 4Ji and J 2 = 2J 2 then 
Ji and J 2 are mapped precisely onto the brackets 7Ti and 112 under the Flaschka transformation. 
It is easy to see that we have a bi-Hamiltonian pair. We define 

/ii = -2(pi +p 2 H hpjv), 

and /12 to be the Hamiltonian 

1 iV-1 
j=l i=l 

Under Flaschka's transformation ([7]), /ii is mapped onto 4(&i + 6 2 + • • • + ^jv) = 4 trL = 4i?i 
and h 2 is mapped onto 2 trL 2 = 4if 2 - Using the relationship 

tt 2 V i/i = TTiV H 2 , 

which follows from part (w) of Theorem [Tj, we obtain, after multiplication by 4, the following 
pair: 

JiV /i 2 = J 2 V /ii. 
The Lenard relations for the eigenvalues translate into 

TTjVXi = AiTTy-iVAi. (12) 

We define the recursion operator as follows: 

ft = JaJT 1 . 

This operator raises degrees and we therefore call it the positive Toda operator. In (q,p) 
coordinates, we use the symbol Xi as a shorthand for X^ 1 , the Hamiltonian vector field of hi 
with respect to the symplectic bracket J\. It is generated, as usual, by 

?Ci = 1^} ^ X\ . 

In a similar fashion we obtain the higher order Poisson tensors 
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We then define the conformal symmetry 

TV 
i=l 

It is straightforward to verify that 

£z Ji — —Ji, 
Czo'h = 0. 

In addition, 

Z (hi) = h u 
Z (h 2 ) = 2h 2 . 

Consequently, Z is a conformal symmetry for J 1; J 2 and hi. The constants appearing in 
Theorem H] are A = — 1, // = and v = 1. Therefore, we end-up with the following deformation 
relations: 

[Z h h k \ = (i + k)h i+k 
Lz.Jk = (k-i- 2)J i+k 
[Zi, Z k \ = (k - i)Z i+k . 

We compute the divergence of the master symmetry Zi, as we will use it in the proof of 
Theorem O It was proved in [5] that Z\ = J2iLi A ?zlr> where A« are the eigenvalues of the Jacobi 
matrix L. Using the definition of divergence we have that D{Z X ) = 2^-Li A i = 2 Ei=i b i = 
~Y^i=iPi = Therefore, 

D{Zi) = \hi. 

The following proposition will also come to use in the proof of Theorem [51 
Proposition 4 

X* = 4Xf j+ \ where f = In(Vdetft) for j>l. (13) 

Proof. 

Similar to the proof of Theorem [2j We note that one uses the fact that the eigenvalues of 7Z 
are the squares of the eigenvalues of L, see [HI |9] . 

Theorem 5 For j >l,Yj is a Hamiltonian vector field given as 

y. = (j - i)xf, where f = ln(Vdet K). (14) 



2' 1) ^ ^ 
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We will prove the theorem for j=l, . . . , 4, and then we will use induction for j > 5. First we 
observe that Y\ = 0, since Ji is symplectic. Using the general form of the tensor J 2 we obtain 
the following: Y 2 = — 2 X)i=i • ^ s i m pl e calculation gives that X^ = — 8^2 i=1 -^-. Thus 
Y 2 = \Xl\ = X J f 2 } using Proposition H 

We have that J 3 = —[Zi, J 2 ] and D(Zi) = Therefore, 



Y 3 = D(J 3 ) = -D[J 2 ,Z 1 ] 

= — ([£)( J 2 ), Zi] — [J 2 , D{Z\)]) 

= -[Y 2 ,Z 1 ] + ±[J 2 ,h 1 ] 



■{xt\z 1 } + -xf 



1 

Using the super-Jacobi identity for the Schouten bracket, the first term equals 

-[J 2 ,[Zi, /]-[/, [J 2 ,^]]. 

Therefore 

= 2X* 

In the last step of the argument we have used that X^ = ^Xj 3 . The proof of the formula 
Y4 = 3Xj! 4 is identical to the one for Y 3 except that we employ the relation [Z 2 , J 2 ] = — 2J 4 . An 
inductive argument based on the same technique that we have used for Y 3 , can also be used to 
show that the result of the theorem holds for j > 5. We omit the details. 

In the theorem that follows, we present an iterative formula that produces all members 
of the modular class in terms of the recursion operator TZ and the modular vector field Y 2 = 
-2£f=ii = (-2,...,-2,0,...,0)*. 

Theorem 6 For j > 2, 



Y j+l = jTV- l Y 2 . 



Proof. 



Y j+1 = jX J / +1 



4/11 4 /ij 4 3 

3 -n j - 1 x l 

3 -n j - l x^ 
4 ™i 

-W~ l AXf 2 
4 ; 

jW~ x Y 2 . 
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5 Summary 



In this paper we study the hierarchy of modular vector fields associated to the infinite family 
of Poisson tensors for the classical Toda lattice equations. We present analytical expressions 
of the modular vector fields both in Flaschka coordinates (a, b) G M 2 ^ -1 , as well as in natural 
coordinates (q,p) G M? N . In both cases, all the members of the infinite modular family, denoted 
by Yj, j > 1, are Hamiltonian. In (a, b) — variables we have that 

Y 3 = ^ln(ai-o A r_i)+(j-l)ln(det(L)) 

where L is the Jacobi matrix of the Lax pair. In natural coordinates (q,p), the modular vector 
field takes the form 

Yj = (j ~ ^^(VtetTi) 

where 1Z is the recursion operator. It is not difficult to show that the term X^ det ^ is the 
projection of the vector field xf 3 . _ — — N under the Flaschka map. The term X? 3 , , makes 

^ J ln(Vdet^) 1 ln(ai- ajv-i) 

its appearance due to the change of coordinates. 

The following properties are proved regarding the behavior of modular vector fields when 
they are applied on the constants of motion Hj, and the Lie derivative of the Poisson tensor nj 
in the direction of the modular vector field Y,- L : 

i) Yi(Hj) = Yj(Hi) 

U) Ly/Kj = -Ly.'Ki. 

In (q,p)— variables, we presented a formula that iteratively produces all members of the 
infinite family of modular vector fields, in terms of the recursion operator. Namely, 

Y 1 = 0, 

N d 

Y j+1 ^jTl j - 1 Y 2 for j>2. 

We conclude with an alternate bi-Hamiltonian formulation of the Toda lattice given by the 
relation X^ = ^(^(l)) > ^ n addition to the existing one Xj/^ 1 = X^ 2 . 
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